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Congestion management is one of the important issues in the deregulated 
power systems. There are several methods to eliminate congestion. Utilizing 
FACTS devices is an appropriate option for large-scale and quick control of 
flows of transmission lines. FACTS devices such as Thyristor Controlled 
Series Capacitor (TCSC) can help to mitigate the transmitting flow of power 


in the congested lines, which leads to an increase in the network loading 





ability as well as reduction of both losses and production costs. Due to 
Keywords: the considerably high price of FACTS devices, it is important to determine 
their optimum location on the network. Accordingly, in this paper, 
the Antlion optimization algorithm (ALO) has been employed to conduct 
a congestion management analysis to determine the optimal location for 
the installation of TCSC, which is simulated on an IEEE 14-bus test system 
subject to satisfy the constraints of the market environment. 
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1. INTRODUCTION 

In recent decades, an important structural reform has been made in many power systems, which has 
changed the power industry from a traditional structure to a restructured modern one. This fundamental 
change in structure and operational rules have become pervasive very soon throughout the world. The former 
is called restructuring of power systems, and the latter is called deregulation. In this regard, the generation, 
transmission, distribution segments and energy services were separated from each other in the first step. 
Then the generation and distribution sectors were divided into several independent companies which may 
have governmental or non-governmental ownerships or may be private equities. 

Subsequently, each one of the generating and distribution companies was allowed to compete with 
other companies in the wholesale electricity market to exchange electrical energy as a seller or a buyer. 
Therefore, the increase in competitiveness of electricity commerce has caused the fair price of electricity 
which is determined based on the supply and demand trade-off mechanism which provides both sides of 
the trade with the level of satisfaction. The reduction of generation costs, the improvement in ancillary 
services quality, and improvement in demand-side satisfaction are other benefits of restructuring in power 
industry. The transmission network is a major obstacle for the deregulation of the power systems because of 
two reasons. The first reason is respect to the technical issues, which implies that it is not possible to separate 
the transmission network like generation or distribution sectors to make it competitive. In addition, 
the requisite of the existence of a proper competition between power providers in supplying electricity is 
the fair and not controlled interconnections across the power grid [1-2]. Although the concept of transmission 
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network congestion already exists in traditional power systems, the term “congestion” has been raised as 
the deregulation of the electricity industry has been started. The meaning of congestion is the use of 
transmission network beyond the permissible operating range. The congestion of transmission lines 
can prevent suppliers from making a new contract, and cause the impossibility of execution of existing 
contracts, as well as curtailments, monopoly of prices in some areas, damage to electrical equipment in the 
system due to unplanned load shedding, price spike, increase in the price of electrical energy in some 
locations, etc. [3-4]. 

There are many ways to reduce transmission congestion. The utilization of FACTS devices is one of 
the most effective ways. Utilizing FACTS devices for congestion management purposes is very useful 
because the limitation on optimum power flow due to the power transmission constraint is basically 
removable by the power flow control. TCSC is one of the most utilized kinds of the FACTS devices that can 
be used to absorb or generate reactive power. TCSC can control the transmission power of the line through 
affecting on the impedance of the target line. An advantage of using this equipment is its quick installation 
compared to the construction of a new transmission line. Therefore, the utilization of these devices to 
eliminate or reduce congestion in the short term is justifiable and sensible. 

In recent years, the determination of the size and the optimal location of these devices in 
the networks have drawn a particular attention to this subject as an optimization problem. Various methods 
have been proposed for finding the optimal size and location of TCSC regard to its generation capacity, 
losses, and costs, and various articles have been published in this context. In [5], a sensitivity analysis of 
the active power flow performance index (PI index) has been used for TCSC and TCPAR optimal locating. 
According to this method, TCSC should be installed in a line that has the most negative sensitivity factor, and 
TCPAR should be installed in a line that has the largest value of sensitivity factor so that the installation of 
FACTS devices in the target line must provide the lowest cost and eliminate the congestion. In [6], 
an approach is introduced to find the optimal location of TCSC subject to reduce the congestion cost (CC) in 
a competitive electricity market considering shadow prices. In this paper, the performance index for TCSC 
placement is a combination of lines' power sensitivity factor and shadow prices. In [7], two options of load 
shedding and utilizing TCSC have been evaluated to manage congestion in a bilateral based power market. 

In [8], the effects of the TCSC on the congestion and prices in an electricity market including 
bilateral contracts are investigated and a LMP-based approach is used. In [9], a study on the optimal location 
of the TCSC for congestion management in an electricity market based on sensitivity analysis and 
considering two goals of reducing the total reactive losses of the system and reducing the active power flow 
performance index (PI index) has been conducted. In [10], the authors have proposed a new method based on 
the total FACTS annual income and cost pertaining to TCSC subject to determine its optimal location in 
order to manage congestion in the restructured electricity markets. In [11], TCSC is utilized in the electricity 
market to improve the ability of the system to transmit more power. In this paper, the sensitivity analysis is 
used for TCSC placement. In [12], multi-objective particle swarm optimization algorithm (MOPSO) and 
sequential quadratic programming (SQP) have been employed with regard to the voltage stability index in 
optimal locating of FACTS devices for congestion management. 

In [13], a new method for the TCSC locating is presented in order to absorb the maximum 
transmittable power by the loads through the network’s branches. In [14], the authors have focused on 
the transmission cost and improving it by installing the TCSC. In [15-16], an optimization method is 
proposed to find the best location to install TCSC regard to maximizing the loadability of the Malaysian 
distribution network based on evolutionary optimization technique. Besides, the increase in the network 
loading with respect to the implications of installing a series of capacitors is investigated using the particle 
swarm algorithm based on birds’ flock behavior. In [17], the particle swarm optimization algorithm (PSO) 
has been used to find the optimal value and the optimal location of TCSC and SVC in order to increase 
the reliability of the system. In [18], the authors have used the bacterial foraging algorithm to optimize 
FACTS devices. In recent years, new evolutionary algorithms such as bat algorithm (BA) [19], glowworm 
swarm optimization algorithm (GSO) [20], gravity search algorithm (GSA) [21], gray wolf algorithm (GWO) 
[22], Shuffled frog leaping algorithm (SFLA) [23], biogeography-based optimization (BBO) algorithm [24], 
and big bang-big crunch (BBBC) optimization algorithm [25] have been widely used to solve optimization 
problems in power system operation and market analysis. In the present article, a new method is proposed 
which is obtained by merging the ALO algorithm and optimal power flow, and this approach is employed to 
determine the optimal TCSC location. The simulation is executed on IEEE 14-bus test system shows ability 
and effectiveness of the proposed approach. 
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2. MODELING AND FORMULATION OF TCSC IN OPTIMAL LOAD FLOW EQUATIONS 
2.1. Static modeling of TCSC 

Figure 1 shows the z model of a transmission line that is installed between the bus i and the bus /. 
Assume that the complex voltage at the i” bus and j™ bus are defined as Vô, and V,Zô, respectively. 


The active and reactive power transmission from bus i to j is represented in (1) and (2): 


P, =V,G,— V, (G, cosô; + B, sin 5, ) 


(1) 


= 2 
Q, =-V?(B, +B 


sh 


)- VV, (G, sin Ô; +B, cos ô) a 


Similarly, the active and reactive power transmission from bus j to i bus is shown as (3) and (4): 


P, =V?G, -V V, (G; cos6, + B, sind, ) 


ji P Sagi i 


(3) 


2 é 
Q; =-V; (B; + Ba) -VV, (G, sin ô; + B; cos ô) (4) 


The transmission line model with incorporating a TCSC that is located between the buses i and j is 
depicted in Figure 2. In the steady state, TCSC is contemplated as a static reactance with the value of -jx-. 
The active and reactive power transmission from the i” bus to the j” bus (and also contrariwise) regard to 
present of TCSC is modelled as (5) to (8): 


P® =V?G, -VV (G; cosd, +B, sinô) 


(5) 
Q; =V; (B; +B„,)-VV, (G, sind, +B, cos 6, ) 6) 
c ay ' Hoe 
P; =V,G, -VV, (G, cos, +B, sin d,) 7) 
OS =-V?(B, +B) +VV,(G, sind, +B, cos6,) (8) 
Bus-i Bus-j Bus-i Bus-j 





Figure 1. The transmission line model Figure 2. The transmission line model with presence 
of TCSC 


Where Gi‘; and Bij are (9-10): 


(9) 





(10) 


The variation of the line’s flow due to series capacitance can be represented as a line without 
a series capacitance with flowing power at the injecting and receiving terminals of the line, as shown in 
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Figure 3. The injected active and reactive power to the bus i (P;, Q;© and to the j™ bus (P;, Qj®) can be 
obtained by (11) to (14) [26-28]: 





PS =V, AG, -V,V; (AG, cosd, + AB, sind, ) aD 
Po =V;AG, - VV, (AG, cosd, + AB, sind, ) ad 
Of =-V?AB, -VV,(AG, cos 6, - AB, sin 5, ) as 
OF =-V7AB, + VV,(AG, sind, + AB, cos 5; ) tid) 
Where AG; and AB; are (15-16): 

AG, = (xer; (3e -24 )) (f2 (r + =) a 
AB, =(~xe(1) = txen (0 +28) +(e) )) ie 


Bus-i Bus-j 





ic je 


Figure 3. TCSC injection model 


2.2. Optimal load flow equations 

The objective function may contain economic, security or environmental aspects of power systems 
that must be solved by a proper optimization algorithm. In recent years, regard to the raising of the concept of 
restructuring and deregulation in the electricity industry, the objective function is mostly defined as 
the minimization of the generation cost (economic aspects of the system) as well as maximization of 
the reliability of the system (security of the system). 


F(x) = min [Se (P; ) ain 


Where F(x) is the objective function that must be optimized, x is the state variables, Ng is 
the number of network generators, c(Pgi) is the generation cost of unit i. In general, our goal is to optimize 
the objective function by a suitable solution and satisfying the prevailing constraints of the system (physical 
constraints, which limit the power generation and availability of transmission lines’ capacities, and 
the constraints imposed on electrical devices used in power grids and system operational strategies). 
If the TCSC is located in the line between buses i and j, then the power balance equation for the nodes i and j 
would be expressed as (18-21): 





P(0,V)- P; +P, +P =0 (18) 
O,(9.V)- Oe, + Op, +Q™ =0 (19) 
P. (0,V)- P; +P, +P =0 (20) 
Q,(0,V)-Q, +Q, +O; =0 (21) 
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Where Pg and Qasi are the active and reactive power in the node i, Pg and Qa; are the active and 
reactive power in the node j, Pp; and Qp; are the active and reactive power consumed by the demand in 
the node i, Pp; and Qp; are the active and reactive power of the loads in the node j, P?@° and Q,7°° are 
the active and reactive injection power by TCSC to the node i, P@° and Q;'°° are the active and 
reactive injection power by TCSC to the node j. The constraints of the problem are also expressed in 
the following (22-26): 





|S, (0.V)< 57" oe 
ies eS P3” (23) 
Qg” < Qo, < QG” (24) 
Vy SVS Veo (25) 
Kp Sapo (26) 


The (22) shows the limitation of the apparent power through the line where Sj is the apparent 
passing power through the transmission line between buses i and j, and Sj"“* is its maximum boundary. 
The (23) and (24) explain the active and reactive power generation limitations so that Pg?" and Pg; imply 
on the minimum and maximum active power generation boundaries in bus i, Qc" and Qg/" are 
the minimum and maximum reactive power generation limits in the bus i. The (25) shows the voltage range 
limitation where V” and V;"“* determine the minimum and maximum limits of the permissible voltage range 
in the bus i. The (26) shows the limitation of the TCSC reactance where Xc""" and Xc"™ are the minimum and 
maximum of TCSC reactance [29-30]. The optimal power flow optimization problem is 
formulated as follows: 





(Sele) Sale P+ Pp HPI) +S Aa (0i -0a +n +0") +r, (5, - sy" )+ Sate (Re -R z) 


i= ij=1 


+5 Ta (P, — pr) 4 ` jp ( ph- P,)+ Y Ta (P; — pr) 4 > a" (vee ba: v,) v) ues (v, 7) 


i=l i=l 


(27) 


Where A and u denote the Lagrangian coefficients of the equality and inequality constraints 
respectively, each of which has an economic interpretation. The most important one is Ap, which is 
the instantaneous price or nodal price or locational marginal price (LMP). Accordingly, by considering 
the presence of TCSC in the network, the overall cost function will be made up of two parts: 

— The cost of power generating at the plant. 
— The cost of investment related to the TCSC. 

Thus, the objective function of the system represents the minimization of the cost of generation 
(economic aspects of the system) and the cost of the TCSC installation: 


F(x) = nin ËG (Pa) + Cioe] (28) 


2.3. Calculation of LMP and congestion analysis 
To calculate the congestion cost in each line, the LMP price difference between two buses must be 
multiplied by the flow of power passing through the line. 


MA, =4-4; (29) 


CC; = A4;P; (30) 
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Where 4; denotes the locational marginal price in the i“ bus, 2; shows the locational marginal price 
in the j™ bus, AA, indicates the marginal price difference between buses i and j, Pi represents the flow of 
passing power through the line i-j, and Cj is the congestion cost of the line i-j. Accordingly, the total 
congestion cost can be calculated as in (31) [31-32]: 


N, 
TCC = > AAP. 
2 Y y 3 1) 


3. ANTLION OPTIMIZATION ALGORITHM 

Antlion optimization (ALO) method is a novel nature-inspired algorithm that is introduced by 
Mirjalili in 2015 [33]. The ALO algorithm mimics the hunting mechanism of antlions in nature. Five main 
steps of hunting prey such as the random walk of ants, building traps, entrapment of ants in traps, catching 
preys, and re-building traps are implemented. There are several different species of ants around the world in 
nature. Antlions belong to the Myrmeleontidae family and Neuroptera order (net-winged insects). 
The lifecycle of antlions includes two main phases: larvae and adult. A natural total lifespan can take up to 3 
years, which mostly occurs in larvae to become an adult ant. The antlion larvae period mostly passed on 
walking routes on sand and leaves to find a good place for building traps. During the hunting process, 
an antlion larva digs a cone-shaped pit in soft sand. As illustrated in Figure 4, after digging the trap, 
the larvae hide underneath the bottom of the cone, and waits for the prey (ants and other kinds of insects) to 
be trapped in the pit. 


Opa ANERIAN; 
Per JAW? 
(GAINS BK WS 





Figure 4. Hunting behavior of antlion 


Once the antlion realizes that a prey is in the trap, it tries to catch it. However, if the prey tries to 
escape from the trap, the antlion intelligently throws sands towards to edge of the pit to slide the prey into 
the bottom of the pit. The mathematical model of ants and antlions is discussed in the following part: 


3.1. Random walks of ants 
The random walks of ants for searching food in the nature can be expressed by (32): 


X(t) = [0, cumsum(2r(t,)-1), cumsum(2r(z, )-1), .... cumsum(2r(t, )-1)] 


(32) 


Where cum sum is the cumulative sum and N is the maximum number of iteration, t shows the step 
of random walk and r(f) is a stochastic function defined as (33). 


D- 1 if rand>0.5 
O= if rand <0.5 (33) 


Where t shows the step of random walk and rand is a random number generated in the interval of 
[0, 1]. The position of ant is presented in the following matrix: 
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A A2 Ang 
M. = Ai Azz a Asa 
Ant T . . . . 
: : ae (34) 
Au A,2 pi Ava 


Where Man is the matrix for saving the position of each ant, Ai; shows the value of the j-th variable 
of i-th ant, n is the number of ants and d is the number of variables. The objective function of each ant is 
saved in Moa matrix. 


F(A Ajas os Ava) 
2 f(A, Aggy ves Aa |) 


M 
OA (35) 
FUA. Aro o Ana ]) 
where f denotes the objective function. 
AL, AL AL, 4 
_ AL, AL,, dia AL, 4 
Antlion : : Py : (36) 
AL, AL AL 


‘n,1 ‘n,2 


Where Maniion is the matrix for saving the position of each antlion, AL; shows the j-th value of i-th 
antlion, n is the number of antlions, and d is the number of variables. The objective function of each ant is 
saved in Mog matrix. Similarly, the objective function of each antlion is saved in Moaz matrix. 


f([ Abs Aline ALa |) 


F([AL» Alby + AL, ]) 
(37) 


f (LAL, ALa AL,a]) 


Moa, = 


In order to keep the random walks inside the search space, a normalizer function is employed 
shown in (38). 


y Biala) 
Eo (38) 





Where a; is the minimum of random walk of i-th variable, b; is the maximum of random walk in i-th 
variable, cf is the minimum of i-th variable at t-th iteration, and d'; indicates the maximum of i-th 
variable at t-th iteration. 


3.2. Trapping in antlion’s pit 
The mathematical model of the trapped ants in the antlion's traps is presented by (39) and (40). 


c= Antlioni, +c! 


(39) 


t > t t 
d; = Antlion; +d (40) 


Where c is the minimum of all variables at t-th iteration, d' indicates the vector including 
the maximum of all variables at t-th iteration, c} is the minimum of all variables for i-th ant, d'; is 
the maximum of all variables for i-th ant, and Antlion‘; shows the selected position of the j-th antlion 
at t-th iteration. 
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3.3. Building trap 
In order to model the antlions’s hunting capability during the optimization process, a roulette wheel 
is employed. This mechanism gives high chances to the fitter antlions for catching ants. 


3.4. Sliding ants toward antlion 

With the mechanisms proposed so far, antlions are able to build traps proportional to their fitness 
and ants are required to move randomly. Moreover, antlions shoot sands outwards the center of the pit to 
slides down the trapped ant that is trying to escape. Equations (41) and (42) expressed the mathematical 
model of ants' slide down and trapped in the antlion's trap. 


d= / I (41) 
d' =d'/I (42) 
where J denote the calculated ratio, obtained by (43) [33]. 

I =10".(t/T) (43) 


Where ¢ is the current iteration, T is the maximum number of iterations, and w is a constant defined 
based on the current iteration and is obtained by (44). 


2 if OIT 
3 if OST 
w=44 if t>0.75T 
: (44) 
5 if 1>0.9T 
6 if t>0.95T 


3.5. Catching prey and re-building the pit 

The final stage of hunt is when an ant (prey) reaches the bottom of the pit and is caught in 
the antlion’s jaw. After this stage, the antlion pulls the prey inside the sand and consumes its body. 
For mimicking this process, it is assumed that catching prey occur when ants becomes fitter (goes inside 
sand) than its corresponding antlion. An antlion is then required to update its position to the latest position of 
the hunted ant to enhance its chance of catching new prey. Equation (45) is proposed in this regard. 


Antlion, = Ant; iff (Ant;)>f (45) 


Where t shows the current iteration, Antlion'; shows the position of selected j-th antlion at t-th 
iteration, and Ant’; indicates the position of i-th ant at t-th iteration. 


3.6. Elitism 

Elitism is an important characteristic of evolutionary algorithms that allows the optimization 
algorithm to select and use the best solution obtained at any stage of optimization process. Since, in 
optimization process antlion is considered as elite, it should be able to affect the movements of all the ants 
(preys) during iterations. Therefore, it is assumed that every ant randomly walks around a selected antlion by 
the roulette wheel and the elite simultaneously. The mathematical model of this behavior is as (46). 


Anti =(Ri, + Re) /2 us 


Where R'4 is the random walk around the antlion selected by the roulette wheel at t-th iteration, R'z 


is the random walk around the elite at t-th iteration, and Ant'; indicates the position of i-th ant 
at t-th iteration [34, 35]. 
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4. SIMULATION AND RESULTS 

In this paper, IEFE 14-bus test system is used for simulation. This system consists of 5 generators, 
11 loads, 17 transmission lines and 3 lines of transformers. The single-line diagram of this network is shown 
in Figure 5. The information on this network is also given in Tables 1 to 3. 

In this method, two cases are investigated: 
— TCSC is not installed in the system and the congestion constraint for transmission lines is imposed. 
— The lines have congestion constraints and TCSC has been added to the system. 

In the first case, in which the system does not have TCSC, the total cost of generation is 
8,156.81 ($/h), and the total congestion cost is 1181.32 ($/h). Thus, the total cost will be 9338.13 ($/h). In the 
second case, after the optimization is executed, the optimal location and size of the TCSC is investigated. 
The line 2 is determined as the best place to install, and X7csc is equal to -10.58. In addition, the total cost of 
generation in the second case is 8061.46 ($/h), the total congestion cost is 816.32 ($/h). Thus, the total cost 
will be 8877, 78 ($/h). Therefore, with the installation of TCSC in the system, the profit will be equal to 
460.35 ($/h). Figure 6 shows the LMP for all buses in the two studied cases. Accordingly, in all buses, aside 
from buses 3 and 8, the LMP price is reduced. In Figure 7, the transmission power through the network lines 
is shown in two cases. 


Table 1. Cost function coefficients and energy sale biddings 








Pus Piz Price ($/MWh) Y B a Bus No. 
332 0 35 0.04743 20 0 1 
140 0 36 0.2391 20 0 2 
100 0 38 0.037 35.4 0 3 
100 0 60 0.02 40 0 6 
100 0 40 0.03 35 0 8 





Table 2. Characteristics of a transmission line 
Line No. BuSsending __ BUSpeceiving Resistance (p.u) _ Reactance (p.u) 








1 1 2 0.01938 0.05917 
2 1 5 0.05403 0.22304 
3 2 3 0.04699 0.19797 
4 2 4 0.05811 0.17632 
5 2 5 0.05695 0.17388 
6 3 + 0.06701 0.17103 
7 4 5 0.01335 0.04211 
8 4 7 0.00 0.20912 
9 4 9 0.00 0.55618 
10 5 6 0.00 0.25202 
11 6 11 0.09498 0.1989 
12 6 12 0.12291 0.25581 
13 6 13 0.06615 0.13027 
14 7 8 0.00 0.17615 
15 7 9 0.00 0.11001 
16 9 10 0.03181 0.08450 
17 9 14 0.12711 0.27038 
18 10 11 0.08205 0.19207 
19 12 13 0.22092 0.19988 
20 13 14 0.17093 0.34802 





Table 3. Data of buses [36] 
Bus No. _ Poen(p.u) — Qoen(p.u) _ Pecon (P-U) Qn (p.u) Bus type — Qmax.gen (P.U)  Qmin.gen (P.u) 








1 2.32 0.00 0.00 0.00 PV 10.0 -10.0 
2 0.4 -0.424 0.2170 0.1270 Swing 0.5 -0.4 
3 0.00 0.00 0.9420 0.1900 PV 0.4 0.00 
4 0.00 0.00 0.4780 0.00 PQ 0.00 0.00 
5 0.00 0.00 0.0760 0.0160 PQ 0.00 0.00 
6 0.00 0.00 0.1120 0.0750 PV 0.24 -0.06 
J 0.00 0.00 0.00 0.00 PQ 0.00 0.00 
8 0.00 0.00 0.00 0.00 PV 0.24 -0.06 
9 0.00 0.00 0.2950 0.1660 PQ 0.00 0.00 
10 0.00 0.00 0.0900 0.0580 PQ 0.00 0.00 
11 0.00 0.00 0.0350 0.0180 PQ 0.00 0.00 
12 0.00 0.00 0.0610 0.0160 PQ 0.00 0.00 
13 0.00 0.00 0.1350 0.0580 PQ 0.00 0.00 
14 0.00 0.00 0.1490 0.0500 PQ 0.00 0.00 
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without TCSC in the IEEE 14-bus test network. without TCSC in the IEEE 14-bus test network 


5. CONCLUSION 

Series FACTS devices can help to increase power system security by controlling the power flow 
passing through transmission lines. Nevertheless, their considerably high cost of them necessitates 
the accurate placement and sizing of them. In this study, a combined method is employed to optimize 
the TCSC parameters and determine its appropriate location to reduce the cost of congestion and to diminish 
the generation cost in power grid. The result of using ALO algorithm and optimal power flow on an IEEE 
14-bus test system is compared in two cases of with and without TCSC. The results show that respect to 
the optimal installation of TCSC (line 2), the total cost will decrease from 9338.13 ($/h) to 8877.78 ($/h). 
This means that the profit will be equal to 460.35 ($/h). Therefore, in congestion management method, 
despite the fact that FACTS devices are expensive, the optimal use of these elements (regard to the best 
FACTS type selection and choosing the best installation location) mitigates the congestion and make a proper 
congestion management possible. Therefore, utilizing these devices for congestion management has a higher 
priority than other methods. 
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